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We identify the factors responsible for the appearance of energy-density inhomogeneities in a 
self-gravitating fluid, and describe the evolution of those factors from an initially homogeneous 
distribution. It is shown that a specific combination of the Weyl tensor and/or local anisotropy of 
pressure and/or dissipative fluxes entails the formation of energy-density inhomogeneities. Different 
cases are analyzed in detail and in the particular case of dissipative fluids, the role of relaxational 
processes as well as non-local effects are brought out. 
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I. INTRODUCTION 

The role of energy-density inhomogeneity in the col- 
lapse of self-gravitating fluids is quite relevant and has 
been extensively discussed in the literature (see 
and references therein). 

Furthermore the Penrose's proposal [9] to define a 
gravitational arrow of time in terms of the Weyl tensor, 
is based on the particularly simple relation between the 
Weyl tensor and density inhomogeneity, for perfect flu- 
ids. However the fact that such a relationship is no longer 
valid in the presence of local anisotropy of the pressure 
and/or dissipative processes, and/or electric charge dis- 
tribution [10l4l2| . explains its failure in scenarios where 
the above-mentioned factors are present (l3l - [l5| . There- 
fore, since the rationale behind Penrose's idea is that 
tidal forces tend to make the gravitating fluid more inho- 
mogeneous as the evolution proceeds, thereby indicating 
the sense of time, it should be clear that all factors aso- 
ciated to energy-density inhomogenity (and not only the 
Weyl tensor) should be present in any definition of the 
gravitational arrow of time a la Penrose. 

However, in spite of its relevance, it is unclear yet how 
different physical phenomena affect and, more precisely, 
produce energy- density inhomogeneities. The pertinence 
of such a question is supported by the arguments above. 

It is the purpose of this work to answer to the above 
mentioned question. It entails in fact three different (but 
related) questions, namely: 

• What aspects of the fluid distribution are related 
(and how) to the existence of energy-density inho- 
mogeneities ? 

• How those factors evolve, starting with an initially 
homogeneous distribution? 
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this last question in turn leads to a third relevant ques- 
tion: 



• Under which conditions an initially homogeneous 
configuration remains homogeneous all along its 
evolution (stability problem)? 



Our procedure will be heavily rely on two differential 
equations relating the Weyl tensor to different physical 
variables. On of them is an evolution equation contain- 
ing time derivatives of those variables, whereas the other 
is a constraint equation implying spatial derivatives. As 
far as we know these equations were first derived by Ellis 
[ill Oj} for configurations without any specific symme- 
try, afterwards they have been reobtained and used by 
different authors (see for example [l(|, fl8j-[2^|)- 

We shall consider general fluid distributions endowed 
with anisotropic pressure and dissipating energy during 
its evolution. The specific physical (microscopic) phe- 
nomena behind these fluid characteristics will be not 
discussed here, instead we shall be concerned only by 
the macroscopic (hydrodynamic) manifestations of those 
phenomena (arguments to justify such kind of fluid dis- 
tributions may be found in [l^|-[25] ). 

Dissipation processes are usually treated invoking 
two possible (opposite) approximations: diffusion and 
streaming out, so we shall do here. In the case of dif- 
fusion approximation a causal transport equation will be 
used, allowing to bring out the effects of pre-relaxational 
phenomena. 

In the specific case of localized configurations we have 
to assume that our fluid distribution is bounded by a 
spherical surface. In order to avoid thin shells on such 
a boundary surface Darmois [26| conditions should be 
imposed. 
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II. ENERGY-MOMENTUM TENSOR, 
RELEVANT VARIABLES AND FIELD 
EQUATIONS 

As is usually assumed in the study of self-gravitating 
compact objects, we shall consider that deviations from 
spherical symmetry are incidental rather than basic fea- 
tures of the process involved. Accordingly we shall re- 
strain ourselves to spherically symmetric fluid distribu- 
tions. 

Thus, we consider a spherically symmetric distribu- 
tion of collapsing fluid, bounded by a spherical surface 
S. The fluid is assumed to be locally anisotropic (prin- 
cipal stresses unequal) and undergoing dissipation in the 
form of heat flow (to model dissipation in the diffusion 
approximation) and null radiation (to model dissipation 
in the free streaming approximation). 

Choosing comoving coordinates inside S, the general 
interior metric can be written 



A dt z + B z dr z + R Z {d9 z + sin^ Odtf ) , ( 1 ) 



where A, B and R are functions of t and r and are 
assumed positive. We number the coordinates x° = t, 
x 1 = r, x 2 = 9 and x 3 = <fi. 

The matter energy-momentum T a p inside £ has the 
form 

T afj = (n + p±)v a v + p ±9a p + (p r - p±)x a x? + q*Vp 

+ V a q, 3 + d a lp, (2) 



where h a p = g a p + V a Vp. 

We do not explicitly add bulk viscosity and/or shear 
viscosity to the system because they can be absorbed 
into the radial and tangential pressures, P r and P±, of 
the collapsing fluid. 

Since we assumed the metric (p} comoving then 

V a =A- 1 5% 1 q a =qB- 1 6f, 
l a = A- 1 S^+B- 1 6^, X a =B- 1 5%, (7) 

where q is a function of t and r satisfying q a = qx a ■ 

From ([5]) with ([7]) we have for the acceleration and its 
scalar a, 



A' ( A' \ 

a 1 = T , o a =o«a ot =^J , 



where a a = ax a , and for the expansion 



1 IB R\ 



(8) 



(9) 



where the prime stands for r differentiation and the dot 
stands for differentiation with respect to t. With ([7]) we 
obtain for the shear (|6]) its non zero components 



oil = 3^ 2(T ' °22 



sin 2 9 3 



R 2 a, (10) 



and its scalar 



where /i is the energy density, P r the radial pressure, 
Pl the tangential pressure, q a the heat flux describing 
dissipation in the diffusion approximation, e the energy 
density of the null fluid describing dissipation in the free 
streaming approximation, V a the four velocity of the 
fluid, x a a un if f° ur vector along the radial direction 
and l a a radial null four vector. These quantities satisfy 



and 



V a V a = -l, V a q a =Q, x a X. a = l, (3) 



X a V a =0, l a V a = -l, l a l a = 0. (4) 



where 



cr H cr a /3 — -a , 



1 B R 



a \b r y 

Then, the shear tensor can be written as 

1, 



OaP = & ( XaX/3 - ^h a/ 3 



(11) 



(12) 



(13) 



The four-acceleration a a and the expansion of the 
fluid are given by 



a a = V aiP V p , e^v a , a . 

and its shear a a p by 



(5) 



= V(a;P) + a{ a Vp) ~ g©^a/3, (6) 



A. The Einstein equations 

Einstein's field equations for the metric ([T]) are given 

by 

G a0 = 8nT a p, (14) 
its non zero components with (p}, ([2]) and ((?]) become 
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SttToo = 8^ + e)A z = [ 2- + - - - ( - 



A 



B 



2— (Ey -o— — -(- 

R + \r) BR \R 



(15) 



8ttT 22 = -^-T 33 = 8irP ± R 2 
sin 9 
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(18) 



B. The mass function 



Let us now introduce the mass function m(t, r) [27 
(see also [28|), denned by 



m- R3 R 23 - R 
m-—R 23 - 



R 



R' , 

B ] ' 



(19) 



Following Misner and Sharp [27|, it is useful to define 
the proper time derivative Dt given by 



and the"i?" derivative Dr, 



1__5 

Adi' 



1 d 



(20) 



(21) 



where R defines the areal radius of a spherical surface 
inside S ( as measured from its area). 

Using (|20|) we can define the velocity U of the collaps- 
ing fluid as the variation of the areal radius with respect 
to proper time, i.e. 



U = D T R. 
Then (fT9l) can be rewritten as 



(22) 



Using (TI5 ]) -([T7 1) with (gDJ) and we obtain from (fT 
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Dttu = —4tt 



and 



Pr --W)U + qE 



U 



D R m = 4ir [jl + q—)R\ 



which implies 



JJ 



m = 4tt / [fi + qj-] R 2 R'dr, 



(24) 



(25) 



(26) 



(assuming a regular centre to the distribution, so m(0) = 
0). Integrating (f2f)|) we find 



with 



fi = H + e, 

Pr =Pr+e, 
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C. Weyl tensor 

The Weyl tensor is denned through the Riemann tensor 
ii^„, the Ricci tensor R a p and the curvature scalar 1Z, 
as: 



Observe that we may also write E a p as: 



(32) 



Finally, using (US]), (HH), (HHJ) with (l9|) and we 
obtain 



1 



1 



The electric part of Weyl tensor is defined by 

with the following non-vanishing components 

2 
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(28) 



(29) 



(30) 



—5-+ ^ + ^ -^2' (3 1 ) 



with 



^= 47 r(A-n)-£, 



n = p r - p_l 



(33) 



III. BIANCHI IDENTITIES AND ELLIS 
EQUATIONS 

As mentioned in the Introduction two differential equa- 
tions for the Weyl tensor will play a central role in 
our work, these two equations originally found by El- 
lis [ljl [13] are here reobtained following the procedure 
adopted in 10]. Beside these two equations we shall also 
need the Bianchi identities, which for the system under 
consideration have two independent components which 
read (see [29| for details): 



(A + Pr) 



B «w~ > i? q'A n JARY „ 

— + 2 (u + P±) — + h 25- = 0, 

B y ^ ' R B q BR 



r 



(34) 



and 



, , A n J B R 



. A' AR' 
» + P ^B +m BR=°- 



(35) 



r 



Finally, the following equations for the Weyl tensor changes in notation) 
may be derived (e.g. see [10] for details, but notice 



[8 - 4tt (A - II)]' = ^ (5 + 4vrF,) + lJhr?^ = ^ [4. + P x ) - £] + 12^, 



(36) 
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[f-47r(A-n)]' = ~(5- 



. „ .BR 3R' 

A,n)-i2, q — = -— 



47r ^ _ ^3 



3m 



- 12 ^AR- 



(37) 



r 



IV. THE TRANSPORT EQUATION 

In the diffusion approximation (e = 0, q = q), we shall 
need a transport equation derived from a causal dissipa- 
tive theory ( e.g. the Miiller-Israel-Stewart second order 
phenomenological theory for dissipative fluids (30l - [33l |). 

Indeed, as it is already well known the Maxwell- Fourier 
law for radiation flux leads to a parabolic equation (dif- 
fusion equation) which predicts propagation of pertur- 
bations with infinite speed (see [34j-[37j] and references 
therein) . This simple fact is at the origin of the patholo- 
gies (38[ found in the approaches of Eckart [39| and Lan- 
dau [40] for relativistic dissipative processes. To over- 
come such difficulties, various relativistic theories with 
non- vanishing relaxation times have been proposed in the 
past [30l - [33l . l4ll |42| . Although the final word on this is- 
sue has not yet been said, the important point is that all 
these theories provide a heat transport equation which is 
not of Maxwell- Fourier type but of Cattaneo type [43| . 
leading thereby to a hyperbolic equation for the propa- 
gation of thermal perturbations. 

A key quantity in these theories is the relaxation time 
t of the corresponding dissipative process. This positive- 
definite quantity has a distinct physical meaning, namely 
the time taken by the system to return spontaneously 
to the steady state (whether of thermodynamic equilib- 
rium or not) after it has been suddenly removed from 
it. Therefore, when studying transient regimes, i.e., the 
evolution from a steady-state situation to a new one, r 
cannot be neglected. In fact, leaving aside that parabolic 
theories are necessarily non-causal, it is obvious that 
whenever the time scale of the problem under consid- 
eration becomes of the order of (or smaller) than the 
relaxation time, the latter cannot be ignored. Indeed, ne- 
glecting the relaxation time ammounts -in this situation- 
to disregarding the whole problem under consideration. 

It is worth mentioning that sometimes in the past it 
has been argued that dissipative processes with relax- 
ation times comparable to the characteristic time of the 
system are out of the hydrodynamic regime. However, 
the concept of hydrodynamic regime involves the ratio 
between the mean free path of fluid particles and the 
characteristic length of the system. When this ratio is 
lower that unity, the fluid is within the hydrodynamic 
regime. When it is larger than unity, the regime be- 
comes Knudsen's. In the latter case the fluid is no longer 
a continuum and even hyperbolic theories cease to be 
realiable. 

Therefore that argument can be valid only if the parti- 
cles making up the fluid are the same ones that transport 
the heat. However, this is never the case. Specifically, 



for a neutron star, t is of the order of the scattering time 
between electrons (which carry the heat) but this fact 
is not an obstacle (no matter how large the mean free 
path of these electrons may be) to consider the neutron 
star as formed by a Fermi fluid of degenerate neutrons. 
The same is true for the second sound in superfluid He- 
lium and solids, and for almost any ordinary fluid. In 
brief, the hydrodynamic regime refers to fluid particles 
that not necessarily (and as a matter of fact, almost 
never) transport the heat. Therefore large relaxation 
times (large mean free paths of particles involved in heat 
transport) does not imply a departure from the hydro- 
dynamic regime (this fact has been streseed before [44j, 
but it is usually overlooked). 

The corresponding transport equation for the heat flux 
reads 



Th ap V^qp. n +q° 



1 



-Kh ap (T tP +Ta p )~KT' 



kT 2 



i|8 

where n denotes the thermal conductivity, and T and 
r denote temperature and relaxation time respectively. 
Observe that, due to the symmetry of the problem, equa- 
tion (|38p only has one independent component, which 
may be written as 

rq = -\nqT 2 (-^) - \rq@A - g(TA)' - qA. (39) 

In the case r = we recover the Eckart -Landau equation. 
In the truncated version of the theory, the last term in 
is absent (see for example 45|), and (f3T)f becomes 



rq + qA=--{TA)'. 



(40) 



V. CAUSES OF INHOMOGENIZATION AND 
STABILITY OF HOMOGENEOUS 
ENERGY-DENSITY CONDITION 

We shall now use the two equations for the Weyl ten- 
sor introduced in section III in order to answer to the 
main questions raised in this work. We shall proceed 
by considering particular cases, in an increasing order of 
complexity 



A. Non— dissipative dust 

In this case P r = P± = q = and from the fact that the 
fluid is geodesic as follows from we can put A = 1. 
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Then equations (1551) and (157)) become 

{£ - 47r/i)' +^§(£- 4^) = 0, 



and 



(41) 



(42) 



From this last equation it follows at once that £ = ■<-)• 
/i' = 0a well known result implying that energy-density 
inhomogeneities in this particular case are controlled by 
the Weyl tensor alone. 

Next, using (fT2"|) and in (T4Tj) , this last equation 
takes the form 



3£i? 



-47T/X(J, 



whose solution may be written as 



47r J fiaR 3 dt 



(43) 



(44) 



where we have chosen the integration function such that 
£(0,r) = 0. 

It is worth noticing that in this particular case (non- 
dissipative dust), conformal flatness and shear-free con- 
ditions imply each other as it follows from an d 
the evolution equation for the shear, which reads (see 

mil) 



17 

¥ 



29er 



-£. 



(45) 



Thus, an initially homogeneous distribution of non- 
dissipative dust will remain so iff it remains conformally 
flat, which in turn implies that it remains shear-free. In 
other words the problem of the stability of the energy- 
density homogenity reduces to the stability of conformal 
flatness or to the stability of the shear-free condition. 



B. Non— dissipative locally isotropic fluid 

The next case in increasing order of complexity, cor- 
responds to a non-dissipative isotropic (in the pressure) 
fluid. Thus we have II = q = 0, P r = P± = P. Then 
equations (IMl) and (I57|) become 



other words, the problem of the stability of the energy- 
density homogeneity reduces to that of the stability of 
the conformal flatness. It is also worth mentioning that 
in this case the vanishing of Weyl implies shear-free, how- 
ever the inverse is no true, since the evolution equation 
for the shear now reads 



a a 
~B ~A 



Next, using using (fT2j) and (J34J) in (|46[) . this last equa- 
tion becomes 



£ 



3£R 



= -47r(/i + P)A<7. 



In the general shearing case we have from (p[9"f 



£ 



4tt / V + P)A<jR?dt 



R? 



(49) 



(50) 



where again we have chosen the integration function such 
that £(0, r) = 0. Obviously, in the case of dust (|50|) 
becomes (|1H) . 

Thus the stability of the energy-density homogeneity 
condition is as in the previos case controlled by the shear 
of the fluid. 

However, as mentioned before, even though all spheri- 
cally symmetric conformally flat spacetimes (without dis- 
sipation) and isotropic fluids are shear-free the inverse is 
not true, and we may assume the fluid to be shear-free 
without that implying conformal flatness. In this latter 
case we have from (|49[) 



£ = 



i? 3 - 



(51) 



where f(r) is an arbitrary function of r satisfying /(0) = 
0. 

Then, if we start with an initially homogeneous con- 
figuration at t = 0, implying £(Q,r) = then f(r) = 
producing £ = for any t and accordingly the energy- 
density homogeneous condition will hold all along the 
evolution. 

If instead we assume a small though non-vanishing ini- 
tial Weyl tensor then it will remain small if the fluid ex- 
pands. However in the contraction case it may become 
relevant for sufficiently small R. Thus, in the case of 
contraction the energy-density homogeneity might not 
be stable in the Lyapounov sense for a shear-free perfect 
fluid. 



(£ - 47r/i)' + P)} = 0, (46) 



and 



(47) 



As in the previous case we have from (1471) that £ = «-> 
fi' = 0, meaning that energy-density inhomogeneities are 
controlled by the Weyl tensor alone, in this case too. In 



C. Locally anisotropic non— dissipative fluid 

We shall now bring out the role of the pressure 
anisotropy in the problem under consideration. For that 
purpose we shall assume q = but n ^ 0. Then equa- 
tions (l36l) and (l37l) become 



(£ - Air^i + AttU) 



3_R 



[£-4tt( M + P ± )] =0, (52) 
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and 



D. Dissipative geodesic dust. 



(£ - 47T/X + 47rn)' + — {£ + 47rn) = 0. (53) 
R 

The first remark emerging from (|53|) is that unlike 
the previous cases, now the responsible for the existence 
of energy-density inhomogeneity is not the Weyl tensor 
alone, but the quantity £ + Anil. Indeed, assuming that 
the fluid distribution is regular at the centre, it follows 
at once from (53]) that £ + Airll = O // = 0. 

Next, using (fT2|) and (l34l) in (|52|) , this last equation 
becomes 

(g + 47rny + (3g + ^ n)fl = -4^ + P r )Aa, (54) 
it 

thus we have an evolution equation for the quantity 
responsible for the energy-density inhomogenity. It is 
worth observing that the quan tity —{£ + 47rII) is referred 
to in previous works 21 1, [46| (notice a difference in no- 
tation with respect to the first of these references) as one 
of the structure scalars (Xtf)- It is defined as follows: 
Let us introduce the tensor X a p defined by: 



X aP =* K 1PS VV 5 = -V ai ep K pf}S VW 5 , (55) 



where R* a p 1& = \Vepi&R a p P anc ^ VepyS denotes the Levi- 
Civita tensor. 

Tensor X a p may also be expressed through its trace and 
its trace-free part, as 

X a p = -X T h a p + X TF (x a Xi3 ~ ^h a p). (56) 

Then from field equations, (f51)]) and using ([25]) (|3"Tj) and 
(l32j) . we obtain (assuming vanishing dissipation) 



X 



TF 



£ - 4ttII. 



(57) 



Therefore the evolution equation ([54]) in terms of Xtf 
reads 



Xtf + 3XT / R = -87rn| + 4tt( M + P r )Aa, 



(58) 



whose general solution is 

4tt J" * hull - (/i + P r )ARa 



X 



TF 



R 2 dt 



R :i 



(59) 



Thus our system is initially homogeneous in the energy- 
density, however as time goes on this condition will be 
affected by the anisotropy of the pressure and the shear 
of the fluid according to (|5T))) . 



Finally, in order to bring out the effect of dissipation 
in the formation of energy-density inhomogeneities we 
shall consider the case of dissipative geodesic dust. Ob- 
serve that in the dissipative case the dust condition does 
not imply that the fluid is geodesic as can be easily seen 
from (1331) . therefore the geodesic condition here is not 
redundant and is assumed, both, for sake of simplicity 
and also in order to isolate the effects produced by dis- 
sipative phenomena alone. It is worth mentioning that a 
great deal of work has been done in order to find exact so- 
lutions descri bing radiating geodesic dust configurations 
(e.g. see [46|-|53l| and references therein) 

Thus we have P r = P± = and A = 1, in which case 
our equations (|3T)j) and (|3"T)) become 

(£ 4^) + 3 -§ {£ - 4vr/i) - 1 -^L = 0, (60) 

and 

(£-47r/i) = — —£ — . (61) 

From the equation above it follows that the quantity that 
now determines the existence of energy-density inhomo- 
geneities is: 



12tt /J" qBRR 2 dr 



(62) 



Indeed, from (1611) is not difficult to find that ji' = 0f> 
* = 0. 

Then, an equation for the evolution of VP can be ob- 
tained from ®, using $F2$ and O : 



Q Ana' AnqR' 3*i? 

* 7 = -Anna H , 

i? 3 B BR R 

with n = 12tt J Q r qBRR 2 dr. 
The general solution of (|63|) is 



(63) 



* = Io^-^R 3 -^ + 4 -^)dt ^ (64) 

where the influence of difcrent factors on the evolution 
of * are clearly indicated. 

We can further transform the expression above by 
noticing that in the case under consideration (1351) can 
be formally integrated to obtain 



0(r) 
B 2 R 2 ' 



(65) 



where <f>{r) is an arbitrary function satisfying 0(0) = 0. 

Thus we see from l|64p that starting from an homoge- 
neous configuration the appearance of density inhomo- 
geneities will depend on two different factors; on the one 
hand the shear of the fluid and on the other dissipative 
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terms. In order to isolate further the influence of these 
latter terms, let us consider the shear-free case. 

If the shear is assumed to vanish, then it is a simple 



J 



matter to see that we may put R 
we obtain from (l64l) 



Br, then using (|65l) . 



Jo 47r/ o r 60(r)r(ir 



47r<p(r)' V 3 

— w — 



(pyr)r 
~R? 



dt 



R ? > 



(66) 



from the above it is apparent that dissipation affect \1/ 
through, both, local and non-local terns. 

Finally it is instructive (always in the shear-free case) 
to bring out the role of relaxational effects in the evolu- 
tion of Thus we shall assume that we are in the purely 
diffusion approximation (e = 0) implying q = q; jj, = jj, 
Thus, from (|35|) it follows 

1 = jj-, (67) 

then combining the above equation with (|40[) we obtain 

rT 



or using ((65 



(r) 



R?T' 



(68) 



(69) 



Feeding back (|69|) into (|66|) we can evaluate the influence 
of r on the evolution of $ for any specific configuration. 



VI. CONCLUSIONS 

We have identified the different aspects of the fluid 
distribution which are responsible for the appearance of 
energy-density inhomogeneities and have found the evo- 
lution equation for the variables representing those as- 
pects. 



For dust and isotropic perfect fluids the relevant fac- 
tor is the Weyl tensor, therefore the stability of the ho- 
mogeneous energy-density condition is equivalent to the 
stability of the conformal flatness condition. 

However for locally anisotropic fluids, inhomogeneities 
are also related to the local anisotropy of pressure. 
The specific variable related to the existence of energy- 
density inhomogeneities was found as well as its evolution 
equation. That specific variable was identified as one of 
the structure scalars. 

Dissipation emerges also as a source of energy-density 
inhomogeneities. In order to isolate the effect of dissi- 
pation from others factors we have considered the dust, 
geodesic case, and as a particular subcase we studied fur- 
ther the shear-free case. The non-local contributions of 
dissipation to the formation of inhomogeneities and the 
role of relaxational effects are the most relevant results 
of this part of the study. 

To conclude, the following remark is in order: The fact 
that the different variables considered here affect the for- 
mation of energy-density inhomogeneities might be in- 
tuitively clear, however the specific form in which these 
variables do so, is far from obvious. Indeed, we may con- 
sider a situation in which any of the above mentioned 
variables are present but are such that their combina- 
tions in (j44|) . (|59|) or (|66|) vanish. In this latter case the 
homogeneous energy-density condition will be stable in 
spite of the fact that all (or some ) of the mentioned 
factors are present. 
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